Truncation Strategies for Optimal
نویسنده
چکیده
yearly review of activities and projects CrosSCutS (triannually): newsletter featuring announcements relevant to our users as well as research highlights in the eld of high-performance computing Speedup Journal (biannually): proceedings of the SPEEDUP Workshops on Vector and Parallel Computing, published on behalf of the SPEEDUP Society User's Guide: manual to hardware and software at CSCS/SCSC To receive one or more of these publications, please send your full name and complete address to: Library CSCS/SCSC via Cantonale CH-6928 Manno Switzerland Fax: +41 (91) 610 8282 Abstract. Optimal Krylov subspace methods like GMRES and GCR have to compute an orthogonal basis for the entire Krylov subspace to compute the minimum residual approximation to the solution. Therefore, when the number of iterations becomes large, the amount of work and the storage requirements become excessive. In practice one has to limit the resources. The most obvious ways to do this are to restart GMRES after some number of iterations, and to keep only some number of the most recent vectors in GCR. This may lead to very poor convergence and even stagnation. Therefore, we will describe a method that reveals which sub-spaces of the Krylov space were important for convergence so far, and exactly how important. This information is then used to select which subspace to keep for or-thogonalizing future search directions. Numerical results indicate this to be a very eeective strategy.
منابع مشابه
Optimal truncation in matching markets
Since no stable matching mechanism can induce truth-telling as a dominant strategy for all participants, there is often room in matching markets for strategic misrepresentation (Roth [25]). In this paper we study a natural form of strategic misrepresentation: reporting a truncation of one’s true preference list. Roth and Rothblum [28] prove an important but abstract result: in certain symmetric...
متن کاملOptimal Weighted Recombination
Weighted recombination is a means for improving the local search performance of evolution strategies. It aims to make effective use of the information available, without significantly increasing computational costs per time step. In this paper, the potential speed-up resulting from using rank-based weighted recombination is investigated. Optimal weights are computed for the sphere model, and co...
متن کاملThe evolutionary origins of Lévy walk foraging
We study through a reaction-diffusion algorithm the influence of landscape diversity on the efficiency of search dynamics. Remarkably, the identical optimal search strategy arises in a wide variety of environments, provided the target density is sparse and the searcher's information is restricted to its close vicinity. Our results strongly impact the current debate on the emergentist vs. evolut...
متن کاملStrategies for Model Reduction: Comparing Different Optimal Bases
Several different ways of constructing optimal bases for efficient dynamical modeling are compared: empirical orthogonal functions (EOFs), optimal persistence patterns (OPPs), and principal interaction patterns (PIPs). Past studies on fluid-dynamical topics have pointed out that EOF-based models can have difficulties reproducing behavior dominated by irregular transitions between different dyna...
متن کاملModel reduction techniques with a-posteriori error analysis for linear-quadratic optimal control problems
The main focus of this paper is on an a-posteriori analysis for different model-order strategies applied to optimal control problems governed by linear parabolic partial differential equations. Based on a perturbation method it is deduced how far the suboptimal control, computed on the basis of the reduced-order model, is from the (unknown) exact one. For the model-order reduction, H2,α-norm op...
متن کاملWavelet Representation of the Potential
In this paper we solve the Laplace equation by solving a boundary integral equation in a wavelet basis. In wavelet basis the discretized equations preserve their good conditioning while their sparse structures speed up iterative solvers. We analyze the behavior of GCR algorithms under wavelet transformations and in non-standard form, the optimal choice of Daubechies lters and truncation strateg...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 1996